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1 Introduction 


The family of charm mesons are predicted by heavy quark effective theory [TJ and 
lattice QCD [2]. The IP states have been well measured by the 5-factories and 
LHCb |3. 31 3 [6]. Evidence for higher mass 5(2600) and 5(2760) states has been 
seen [5, 6J. Only natural spin-parity resonances (J F = 0 + , 1 , 2 + ,...) contribute in 
B —» Dt s -)hh' decays where h and h' are kaons and pions. In 2014 LHCb published 
results from a Dalitz plot analysis of 5° — y D°K~ti + decays, which included the first 
observation of the 5 s i(2860)~ and 5 s3 (2860)~ mesons [7, 8]. These states are thought 
to be members of the D s ID family p, [10]. It is therefore interesting to explore 5 
meson spectroscopy to find and identify new states to compare their properties with 
the theory predictions. Three analyses are presented, using B + —> D~K + n + [11] . 
B° —> D°TT + n~ [12] and B° —> D°K + n~ decays [13] . 

2 Dalitz plot analysis of B + —> D~K + n + decays 

The first observation of the decay B + —> D~K + n + , with 5~ —> K + n~'ji ~, is made 
using the topologically similar decay B + —>■ 5~7t + 7t + as a normalisation channel [IT] . 
Event selection is based on a neural network used to reduce combinatorial background. 
Candidates in the signal and normalisation channels are shown in Fig. [lj with fits 
used to extract the signal and normalisation channel yields overlaid. Accounting for 
the selection efficiencies gives the branching fraction ratio 

B(B+ D~K+ 7T+) 

B{B+ D -^ +) = 0 0720 ± 0 0019 ± 0 0021 • t 1 ) 

where the uncertainties are statistical and systematic, respectively. Using the known 
value of B(B + —> D~n + ir + ) = (1.01 ± 0.05) x 10^ 3 [T^ gives 

B(B + ->• D~K + n + ) = (7.31 ± 0.19 ± 0.22 ± 0.39) x 10' 5 , (2) 

where the third uncertainty is from B(B + —y D~h + ti + ). 

The Dalitz plot analysis is performed on candidates in the B mass window 5239.4 - 
5317.1 MeV (natural units are used throughout), with about 2000 signal candidates 
and a purity of approximately 93%. In B + —* 5~A' + 7t + decays, resonances are 
only expected to appear in m(57r), allowing angular moments from the Legendre 
polynomials to be used to guide the amplitude model. The moments study showed 
no evidence of structures above spin 2. The components included in the amplitude 
model are shown in Table [T] 

The amplitude fit is performed with the Laura++ package [15] using the isobar 
formalism nanzmsi, with histograms to describe backgrounds and signal efficiency. 
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Figure 1: Fits to the B candidate invariant mass distribution for (left) B + —> D 7r + 7r + 
and (right) B + —> D~K + tt + candidates. Components are described in the legend. 


Resonance 

Spin 

DP axis 

Model 

Parameters 

£>*(2400)° 

0 

?7T 2 (£>7T ) 

RBW 

m = 2318 ± 29MeV, T = 267 ± 40MeV 

£>2(2460)° 

2 

m 2 (Dn) 

RBW 

Determined from data 

£>}(2760)° 

1 

m 2 (£>7r) 

RBW 

Determined from data 

Nonresonant 

0 

m 2 (£>7r) 

EFF 

Determined from data 

Nonresonant 

1 

m 2 (£>7r) 

EFF 

Determined from data 

£>*(2007)° 

1 

m 2 (Dn) 

RBW 

m = 2006.98 ± 0.15MeV, V = 2.1MeV 

Bf 

1 

m 2 (DK) 

RBW 

m = 5325.2 ± 0.4MeV, Y = O.OMeV 


Table 1: Components of the B + —> D~K + ir + amplitude fit model. RBW and EFF are 
the relativistic Breit-Wigner function and exponential form factor, respectively. Terms with 
subscript v are virtual components, where the resonant pole mass is outside of the Dalitz 
plot boundary. 

For the full fit results see Ref. [nj. Figure [2] shows the fit projection in m(£>7r). Sig¬ 
nificant contributions are seen from the £>q( 2400)°, £>2(2460)° and £>}(2760)° states, 
where the spin of the latter is determined to be 1 for the first time. Other spin 
hypotheses are rejected with high significance (> 6cr). The mass and width for the 
£>^(2760)° and £>2(2460)° resonances are found to be 

?n(£>;(2460)°) = (2464.0 ± 1.4 ± 0.5 ± 0.2) MeV, 

r(£>;(2460)°) = (43.8 ±2.9 ±1.7 ±0.6) MeV, 

m(£>i(2760)°) = (2781 ± 18 ± 11 ± 6) MeV, 

r(£>*(2760)°) = (177 ± 32 ± 20 ± 7) MeV, 

where the uncertainties are statistical, experimental systematic and model dependent 
systematic, respectively. 
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Data D*(2460)° 

Total - D*(2760)° 

Background BJ 0 

D; v (2007)° Nonresonant S-wave 

DJ(2400)° - Nonresonant P-wave 


Figure 2: Projection of the amplitude fit in m(Dir) for B + —> D K + ir + candidates. 
Components are as described in the legend. 

3 Dalitz plot analysis of B° -> decays 

The amplitude analysis of the B° —> D° 7r + 7r~ final state is performed using the 
D° —y K + n~ decay [12]. With larger data samples this channel can be used to measure 
cos(2/3) and sin(2/3) [111], SO], where j3 is an angle of the unitarity triangle. Resonant 
structures are expected in both m(Dn) and m(mi). Combinatorial background is 
removed using a Fisher discriminant multivariate selection. Figure [3] shows the B 
candidate invariant mass distribution of selected candidates. The signal window used 
in the amplitude analysis is 5250-5310 MeV. It contains about 10000 signal candidates 
with a signal purity of around 98 %. 



Figure 3: Fit to the B candidate invariant mass distribution for B° —>• H°7r + 7r decays. The 
data are black points and the fit and backgrounds are shown in red and green, respectively. 

Two amplitude fits are performed, using the isobar model and a K-matrix ap¬ 
proach (2T, 22] for the nn S-wave contribution. The resonances included in the mod¬ 
els are shown in Table [2] Projections of the isobar model fit are shown in Fig [4j 
see Ref. [12] for the K-matrix results. The charm resonances Dq (2400)“, _D2(2460) - 
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and D*j( 2760)“ are found to be significant and the D*j( 2760)' state is determined, 
with high significance, to be spin 3 for the first time. It is interesting to note that in 
the B + —> D~K + tt + analysis the D*j( 2760)° was found to be spin 1. This suggests 
that there could be two overlapping states, as was seen in the D s meson family in 
B J? —y D°K~ 7r + decays [71, Ej- The masses and widths of the charm resonances from 
the isobar model fit are 


m(T>o(2400)“) 

T(D*(2400)-) 

r(D;(2460)-) 
m(D* 3 ( 2760)“) 
r(D*(2760)“) 


(2349 ± 6 ± 1± 4) MeV, 
(217 ± 13 ± 5 ± 12) MeV, 
(2468.6 ± 0.6 ± 0.0 ± 0.3) MeV, 
(47.3 ±1.5 ±0.3 ±0.6) MeV, 
(2798 ± 7 ± 1± 7) MeV, 
(105 ± 18 ± 6 ± 23) MeV, 


where the uncertainties are statistical, experimental systematic and model dependent 
systematic, respectively. Good agreement is seen between the isobar model and K- 
matrix fit results. 


Resonance 

Spin 

Model 

m r (Me 

:V) 


r 0 (MeV) 

D 0 7i~ P-wave 

1 

(Ej 



Floated 

.Dq(2400) - 

0 

RBW 



Floated 

±>2 (2460)“ 

2 

RBW 



Floated 

D}( 2760)- 

3 

RBW 



Floated 

p(770) 

1 

GS 

775.02 

± 

0.35 

» 149.59 ±0.67 

cu(782) 

1 

m 

781.91 

± 

0.24 

8.13 ±0.45 

p(1450) 

1 

GS 

1493 

± 

15 


427 ±31 

p(1700) 

1 

GS 

1861 

± 

17 


316 ±26 

/ 2 ( 1270) 

2 

RBW 

1275.1 

± 

1.2 


185.1 

7T7T S-wave 

0 

K-matrix 




ua 


/o(500) 

0 

na 




m 


/o(980) 

0 

FLT 






/o(2020) 

0 

RBW 

1992 

± 

16 


442 ± 60 

Nonresonant 

0 

na 




m 



Table 2: Components of the isobar and K-matrix amplitude fit models to B° —> H°7r + 7r~ 
decays. GS is the Gounaris-Sakurai function and FLT is the Flatte shape. See Ref. [12] for 
more details. 
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Figure 4: Projection of the isobar model fit in (left) m(Tiir) and (right) m(Dir) for B° —> 
D°Tr + ir~ candidates. The components are (black) data, (blue) isobar fit, (green) 7r + 7r - P- 
and D-wave, (pink) Dtt contributions, (red) 7r + 7r _ S-wave and (grey) background. 

4 Dalitz plot analysis of B° —> D°K + tt~ decays 

An amplitude analysis of B° —>■ D°K + tt~ decays with D° —>• K + tt~ is presented [13]. 
The goal of studying B° —* DK + tt~ decays is to measure the unitarity triangle angle 
7 , as outlined in Refs. [23U2HJ- It can also be used to access the same charm resonances 
as B° —y D 0 7i + 7r~ decays, although the available statistics are smaller. Contributions 
also appear in the m 2 (Kn) axis of the Dalitz plot. 

The event selection is based on a neural network to distinguish between signal 
and combinatorial background. The B candidate mass distribution of selected events 
is shown in Fig. [5j overlaid with the fit used to determine the signal and background 
yields. Events in the signal region, defined as 5248.6-5309.1 MeV, are selected for 
the Dalitz plot fit. There are approximately 2500 signal candidates with a purity of 
around 75 % in this window. 




Figure 5: Fit to the B candidate invariant mass distribution for B° —> D°K + ir candidates 
with (left) linear and (right) log y-axis scales. Components are described in the legend. 
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The amplitude fit contains contributions from the terms shown in Table [3] and is 
performed using the Laura++ package [15] with the isobar formalism. Backgrounds 
and efficiency corrections are both accounted for in the fit. For the full results of 
the amplitude fit see Ref. [13]• The projections of the amplitude fit in m(DTr) and 
m(Kn) are shown in Fig. [6] (left) and (right), respectively. 


Resonance 

Spin 

DP axis 

Model 


Parameters (MeV) 

iC*(892)° 

1 

m 2 (Kn) 

RBW 

m 0 = 

895.81 ±0.19, T 0 = 47.4 ±0.6 

iC*(1410)° 

1 

7n 2 (Kn) 

RBW 

m 0 

= 1414 ± 15, r 0 = 232 ± 21 

K:* (1430)° 

0 

7n 2 (Kiv) 

LASS 


Determined from data 

Ji 2 *(1430)° 

2 

7n 2 (Kiv) 

RBW 

m 0 : 

= 1432.4 ±1.3, r 0 = 109 ±5 

_Dq(2400)~ 

0 

77i 2 (Dtt) 

RBW 


Determined from data 

D*(2460)- 

2 

771 2 (D'k) 

RBW 


Determined from data 

Nonresonant 

0 

tti 2 (Dtt) 

dabba 


Fixed 

Nonresonant 

1 

m 2 (Dir) 

EFF 


Determined from data 


Table 3: Components included in the B° —> D°I\ + tt amplitude fit model. More details 
on the dabba and LASS models can be found in Ref. [ T3l . 



Data - Total fit 

Kji S-wave (LASS) - K*(1430)° 

Du S-wave (EFF) - Dji P-wavc 



. AC*(892)° - fC*(1410)° 

- ' D* (2400) - D 2 * (2460) 

(EFF) D* Background 


Figure 6: Projection of the amplitude fit in (left) m(Dir) and (right) m(A±) for B° —> 
D°K + tt~ candidates. Components are described in the legend. The dip in the right plot is 
due to a veto at m(D°). 

The charm resonance results are in agreement with the B° —> analysis. 

Due to the lower statistics available no contribution is seen at m^Dir) & 2760 MeV. 
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The masses and widths of the states Dq( 2400) and ±>2(2460) are reported to be 


m(L>o(2400)-) 

r( J D*(2400)“) 

771 (^( 2460 )-) 

r(T>;(2460)-) 


(2360 ± 15 ± 12 ± 28) MeV, 
(255 ± 26 ± 20 ± 47) MeV, 
(2465.6 ± 1.8 ± 0.5 ± 1.2) MeV, 
(46.0 ±3.4 ±1.4 ±2.9) MeV, 


where the uncertainties are statistical, experimental systematic and model dependent 
systematic, respectively. These results are in agreement with those from the B° —y 
D°7i + n~ analysis but are less precise. 


5 Summary 

The latest results on charm spectroscopy from Dalitz plot analyses of B meson decays 
at LHCb are presented. First observations are made of the ±>j(2760)° and ±>3(2760)“ 
mesons. Larger data samples are needed to determine whether or not the isospin 
partners of these states can be seen. 
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